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Abstract. We analyse the observational characteristics for a set of Long Duration Events 
(LDE) using Reuven Ramaty High Energy Solar Spectroscopic Imager {RHESSI). Excel- 
lent energy resolution (1 keV) of HXR images reconstructed with RHESSI allowed us to 
perform imaging spectroscopy of coronal loop top sources (LTS) observed for hours after 
maximum of the flare. We found that the sources are large structures. The diameter of 
an LTS is changing with time and energy. High temperature component can be observed 
for few dozen hours after the maximum. For some time intervals clear non-thermal com- 
ponent was observed in spectra, but there was no significant footpoint emission at that 
time. All thermal and non-thermal emissions originate from the LTS. The obtained ob- 
servational char acteristics were used in c alculating the energy balance within LTS, which 
are presented in iKolomahski et al.l l^Ollf lfthis issue). 
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1. Introduction 



Veselovsky Prokhorovl (120081 ) analysed duration times for over 40000 flares 
observed by GOES satellites for years 1976 - 2006. The authors found that 
the histogram of durations can be easily described with the use of lognor- 
mal distribution. With such a distribution it is possible to divide flares into 
three groups: impulsive - shortest duration, typical - located close to the 
maximum of the histogram and long duration events (LDE) - located to the 
right of the maximum. Such a division is arbitrary. There is no difference in 
the physical mechanisms responsible for each group as these flares belong 
to one distribution with no significant fluctuations. As at the basis there 
is no difference between physics of solar flares from the three groups then 
the investigation of LDEs gives better insight into the nature of the energy 
release process with present instrumentation. 
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The decrease of brightness of LDEs observed in soft X-rays (SXR) may 
last from several hours to more than a day. Much insight into the na- 
ture of LDEs was made by ultraviolet and X-ray observations, during the 
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the most important conclusions is that without the continuous energy input 
during the whole decay phase the LDEs would decay much faster than it is 
observed. 

The whole X-ray emission of LDEs is located in loop-top sour ces (LTSs) 



which form before t he flare maximum and last for many hours (e.g. lFeldman et al 



19951 : iKolomahskil . l2007al ). LTSs are very promising features for an analy 



sis of the energy release during the de cay phase as they s hould be close 
to the primary energy r e lease site (e.g. IKopp &; Pneumanl . Il976l : IShibatal . 
19991 : Iffirose et all l200ll : iKarlickv fe Bartal . I2006h 7 Moreover, the presence 



of LTSs during the whole flare decay-phase requires continuous energy re- 
lease and some restriction m echanism effici e ntly p reventing outflow of mass 



and energy from them (see IVorpahl et all 119771 ) . Observations show that 



without such mechanism LTSs would rapidly lose energy by r adiative and 



conductive proces s es and should vanish within minutes (e.g. IJiang et al. 
20061 : IKolomahskil . l2007bh . 

Hard X-ra y LTSs were rarely observed for hours after flare maximum. 



Masuda et al .1 (119981 ) observed HXR sources up to 30 minutes after flare 
maximum. The sources wer e large (1 — 2 arcmin i n dia meter). Two other 
LDEs were investigated by lHarra-Murnion et al.l (jl998| ). The sources ob- 
served grew with time and had diameters ~ 20 — 45 arcsec. In one of 
these flares the HXR source was visible up to 3 hours after flare maxi- 
mum. The analysis of three LDEs with the use of Yohkoh data was made 
recently by Kolomahski (2007a, b). The author reported that HXR sources 
(14 — 23 keV) were observed for about 50 minutes after flare maximum. 
The HXR s ource existin g more than 1 hour after flare maximum was also 
reported by iKhan et al.l (|2006[ ) . In all these cases the HXR emission was 
observed close to and above the SXR emission source. 

First investigation of l ong p ersisting LTS based on RHESSI data was 
made by iGallagher et al.l (120021 ). The HXR emission was observed in the 
12 — 25 keV range some 4 hours after flare maximum and almost 11 hours 
in the 6 — 12 keV range. The LTS was large and its altitude increased at 
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Table I: List of analyzed flares. (1) - Date, (2) - GOES class, (3) - Heliographic coordi- 
nates, (4) - GOES maximum [UT], (5) - GOES decay time [h] . 
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a speed that gradually declined from 10 to 1.7 kms -1 . The higher energy 
emission (12 — 25 keV) was located above the lower energy emission (6 — 
12 keV and 3 — 6 keV) and the whole HXR emission was located above 
tops of the loops observed in the EUV range (TRACE 195 A). Anothe r 
long persisting LTS has been reported recently ( Saint-Hilaire et al. , 20091 ). 
The authors observed an LTS in the energy range 6 — 12 keV even 12 hours 
after the flare maximum. The source observed was spatially correlated with 
the FeXVIII emission observed by Ultra Violet Coronograph Spectrometer 
(UVCS) onboard SOHO. 

Yohkoh/HXT and RHESSI data have provided useful information about 
LDEs, but only in one channel significant emission was observed. Poor en- 
ergy resolution of Yohkoh/HXT was caused by the instrument character- 
istics. Despite the fact that RHESSI has excellent energy resolution, the 
whole previous analysis was restricted to wide energy intervals which still 
do not distinguish the actual nature of long persisting LTS. Here we present 
the investigation made with the use of RHESSI images reconstructed in very 
narrow (1 keV) energy intervals. RHESSI is very useful for analyzing weak 
sources even with 1 keV energy resolution due to its high sensitivity. Thus, 
it gives an opportunity to determine the actual nature of emission of LTSs 
(thermal, non-thermal). Using the images we estimate physical parameters 
of LTSs through imaging spectroscopy. 
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Figure 1: RHESSI light curves of the entire event no. 4. Vertical lines represent borders of 
the spacecraft night (N) and the SAA (S) periods. Arrows mark time intervals in which 
the analysis was performed. 



2. Observations and data analysis 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 
is a rota ting Fouri e r ima ger with nine detectors made of pure germanium 
crystals ( Lin et all . 120021 ). The detectors record energy and time of arrival 
for each detected HXR photon. Pairs of grids placed ahead of the detectors 
and the rotation of the whole satellite (4 s period) cause modulation of 
HXR radiation coming from actual solar sources. Using several available 
algorithms we can reconstruct the spat ial distribution of the HXR emission 
from the time-modulated light curves (IHurford et all |2002[ ) . 

In our analysis we focused on the very late phase of LDEs when the 
emission is extremely weak, so some care is needed in the analysis. It is a 
great difficulty to reconstruct a RHESSI image when the emitting source 
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Figure 2: TRACE images (171 A) with RHESSI contours (6-7 keV) for the 2004 January 
5 flare. 



is weak. Usually RHESSI image reconstruction is performed for detectors 
Nos. 3-6, 8 and 9. Depending on the chosen reconstruction algorithm and 
weights this set of detectors giv es a spatial resolution of the order of 7-9 
arcsec (lAschwanden et all 120021 ), However, we were not able to reconstruct 
any reliable image for the very late phase of an LDE decay. The weak signal 
was not the cause since taking integration times up to 8 minutes we were 
able to collect enough counts. 

The problem was resolved in the following way. First, we reconstructed 
images for single grids and a wide field of view (8 arcsec pixel, image size 
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Figure 3: Spectra obtained for the 2004 January 5 flare. Times are the same as in Fig. [5] 
Thermal component and 6.7 keV line fits are shown. The excess observed for energies 
greater than 10 keV has been fitted with power-law function. Obtained slopes are high, 
of the order of 10. 



256x256 pixels) using the back-projection algorithm (jHurford et all 120021 ) . 
Next, we used these images for determining which grid provides us with 
reliable images of the source. This step is crucial if we remember that when 
a source size is comparable to the resolution of a particular grid then the 



detec tor records very weak or no modulation of the signal (jHurford et al. 



20021 ). In such a case there is only noise in the fine grids even if the count 



rates are high. This caused fail in reconstruction process performed with 
the use of the standard set of grids. 
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Figure 4'- Time evolution of observational characteristics obtained for the 2004 January 
5 flare. 



In the next step we used only selected coarse grid s and reconstructed im- 
ages with PIXON algorithm ( Puetter fe Yahi3 . 1999I . and references within). 
Since integration times were rather long (minutes) we used stacked mod- 
ulations 0. The energy ranges chosen for reconstruction were as narrow as 
possible, i.e. 1 keV hence it was possible to obtain reliable fits to the spec- 
tra derived from the RHESSI images. The image spectroscopy we performed 
has an advantage in comparison to "standard" spectroscopy based on the 
fluxes measured for the whole Sun. Namely, for the rotating Fourier imager 
the background photons do not influence the modulation profile. Thus, the 
analysis of very weak fluxes can be done more precisely especially for mo- 
ments of passage through radiation belts when the background is changing 
rapidly with time. 



1 http: / / sprg.ssl.berkeley.edu/~tohban/nuggets/ ?page=article&article_id=39 
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We did error analysis taking the advantage of slow evolution of LDEs. 
Namely, we divided each analysed time interval into three intervals. Images 
were reconstructed for each interval and the observational characteristics 
were obtained. Next, we calculated a mean value of each observational char- 
acteristic and its standard deviation. We treated the standard deviation as 
the uncertainty of the obtained characteristics. This method was used for 
each parameter: size, location, emission measure and temperature. 

3. Results 

We investigated a set of 9 flares listed in the Table |TJ Fig [T] shows RHESSI 
light curve for energy intervals: 6 — 12, 12 — 25 keV and 7 — 20 MeV for the 
2004 January 5 flare. The last interval was chosen because it shows clearly 
passages through radiation belts. The flare HXR emission in the 6—12 keV 
band was observed for 15 hours after flare maximum recorded at 3:45 UT. 
We marked analysed time intervals with arrows. As mentioned above each 
interval was divided into three subintervals of the same length. 

Reconstructed HXR sources were compared to TRACE images taken 
with the use of 171 A filter. An example of position of HXR sources (6 — 
7 keV) overlaid on TRA CE images is shown in Fig. [2] The HXR sources were 
observed above EUV loops and systematically changed their height with 
velocity of about 2 km/s. We observed no significant footpoint emission in 
the 2004 January 5 flare and in other events also. 

LDEs are characterized by slow decrease of brightness caused mainly by 
slow decrease of temperature. Thus, we tried to define an objective method 
to define the length of the decay phase. For this purpose we fit a Poisson 
distribution to the observed temperature changes: 

T(t) = Ae~r (1) 

where A is a normalization factor and r is a characteristic time. We assumed 
that the asymptote of this distribution is 2 MK. The characteristic time was 
estimated for a group of typical, short duration events and we found that it 
had values of the order of minutes. For the analysed group of LDEs r was 
never less than 0.5 hour (see Table ITT)) . 

Obtained observational characteristics of the 2004 January 5 flare are 
shown in Fig. [4] Both, temperature and emission measure, decrease with 
time. For temperature we observed very fast decrease shortly after flare 
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Table II: Observational characteristics of analysed LDEs. (1) - Date, (2) - Temperature 
[MKj, (3) - Emission measure [10 47 cm" 3 ], (4) - Radius [Mm], (5) - Altitude [Mm], (6) 
- Density [10 9 cm -3 ], (7) - Characteristic time of temperature decay [hj. 
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maximum. Next, when temperature reached value 12 MK, this decrease 
was significantly slower. In the time interval of 6:00 - 18:00 UT temperature 
changed for 3 MK only. 

The evolution of the radius of the observed source shows two phases. At 
first the source expanded from r = 16 arcsec to r = 40 arcsec at 12:00 UT. 
After 15:00 UT we observed that source started to diminish to the radius 
of the order of 28 arcsec. This observation is very important from method- 
ological point of view. Namely we chose grids dynamically depending on 
observed modulation. This result shows that we do not miss information 
about the actual size when using only coarse grids. 

We observed non-thermal, power-law component in fitted spectra (Fig. [3]). 
These have to be treated very carefully since slopes are large, close to the 
value of 10. However, even very steep power-law function is still different 
from thermal one. We observed non-thermal components in six analysed 
flares (among 9). It should be mentioned that non-thermal emission was 
observed in loop-top sources. We did not observe footpoint sources. 
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4. Summary 

We investigated set of 9 LDEs well observed by RHESSI. In Table [H] we 
summarized all obtained observational characteristics. For each parameter 
(except characteristic decay time) we show two values. The first was ob- 
tained shortly after flare maximum and the second was obtained for the 
last time interval. Such intervals show ranges of values that were observed. 
The main conclusions of our work are as follow: 

1. The sensitivity of the RHESSI detectors allows us to observe the decay 
phase in HXR emission up to several hours after flare maximum. 

2. Observed loop-top sources are large structures and they usually grow 
with time. Sometimes we observed decrease of radius after initial 
growth. 

3. HXR emission sources are located above loops observed in the EUV 
range. We did not detect any footpoint emission. 

4. Obtained values give strong constraints for models of solar flares since 
high temperatures demand continuous energy release. 

5. The observed non-thermal component is another proof for continuing 
energy release. However it must be stressed that observed power-law 
components have very steep spectra with indices close to the value of 
10. 
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